Role of T=0 pairing in Gamow-Teller states in N=Z nuclei 
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Gamow-Teller (GT) states in N=Z nuclei with the mass number A from 48 to 64 are studied by using Hartree-Fock-Bogoliubov + 
quasi-particle random phase approximation (HFB+QRPA) with Skyrme interactions. The isoscalar spin-triplet (T=0,S=1) pairing 
interaction is taken into account in QRPA calculations. It is found in the context of SU(4) symmetry in the spin-isospin space that 
the GT strength of lower energy excitations is largely enhanced by the T=0 pairing interaction which works cooperatively with the 
T=l pairing interaction in the ground state. A two-peaked structure observed recently in (p, n) reaction on 56 Ni can be considered 
as a manifestation of the role of T=0 pairing in the GT collective states. 
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The most prominent evidence of pairing correlation in nu- 
clei is found in the odd-even staggering in binding ener- 
gies and the gap in the excitation spectrum of even-even nu- 
clei in contrast to the compressed quasi-particle spectrum in 
odd- A nuclei HI 0, H. There are also dynamical effects of 
pairing correlations seen in the moment of inertia associated 
with nuclear rotation and large amplitude collective motion 
;H H H. The Hartree-Fock (HF)+BCS method and Hartree- 
Fock-Bogoliubov (HFB) method have been commonly used to 
study the ground state properties of superfluid nuclei in a broad 
mass region . For the study of excited spectra, quasi- 

particle random phase approximation (QRPA) has often been 

adopted as a basic method EH Ei El El. 

The strong attraction between nucleons is the basic ingre- 
dient for the pairing correlations. So far, the pairing inter- 
actions of like-nucleons with the isovector spin-singlet (T=l, 
S=0) channel has been mainly discussed. In fact, the attraction 
between protons and neutrons is even stronger in the isoscalar 
spin-triplet (T=0, S=l) channel 114il, which gives rise to the 
deuteron bound state. However the role of T=0 pairing is lim- 
ited in nuclei because of large imbalance between neutron and 
proton numbers, and also the two-body spin-orbit interaction 
which breaks the S=l pair more effectively than the S=0 pair 
J3, E3. EUl- Nevertheless, the isoscalar pairing causes extra 
binding energies in nuclei with N=Z and has been considered 
as one of the origins of the Wigner energy lfl7ll . 

Gamow-Teller (GT) states have been studied both experi- 
mentally and theoretically intensively in the last three decades. 
Many interesting nuclear structure information has been re- 
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vealed by these studies, for example, the quenching of sum 
rule strength 111 811 and the role of GT strength in the astrophysi- 
cal processes such as neutrino-nucleus reactions [19]. Because 
of recent development of modern radioactive beam accelerator, 
it becomes feasible to observe GT states in exotic nuclei near 
the proton and neutron drip lines. Recently, the GT transition 
strength was studied in a N=Z nucleus 56 Ni which has an im- 
portant impact on late stellar evolution through electron capture 
and p decay j20ll . Although the collective GT state is mainly 
built of charge exchange particle-hole excitations, the low-lying 
strength responsible for (3 decay involves neutron-proton {np) 
particle-particle type excitations which could be sensitive to the 
T=0 pairing interaction I10ll2lll . Intuitively, the low-energy GT 
excitations will be enhanced by the following mechanism: the 
T=l pairing correlation provides partially occupied proton and 
the neutron orbitals near the Fermi surface. These orbitals will 
accept the additional proton excitations from the neutron or- 
bits with the same orbital / quantum number. Then, these (np) 
pair interacts through the T=0 pairing interaction and enhances 
GT strength. These physical mechanisms may implement the 
super-allowed GT transition between the same multiplet mem- 
bers of SU(4) symmetry in the spin-isospin space which never 
happened when a neutron orbital is completely full and the cor- 
responding proton one is empty. In this letter, we study GT 
states in several N=Z nuclei focusing on the cooperative role 
of T=l and T -0 p airing interactions in the GT transitions. In 
Refs. B23LI24L 12511 . we applied successfully the self-consistent 
HF+RPA model to calculate spin and isospin dependent exci- 
tations in several closed shell nuclei. As an extension to open 
shell nuclei, the HFB+QRPA method is used to calculate GT 
states with realistic Skyrme interactions. 
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Figure 1: (Color online) Gamow-Teller strength in Ni obtained by 
HFB+QRPA calculations with the Skyrme interaction SGII+Tel. The excita- 
tion energy is referred to the ground state of 56 Ni. The T=0 pairing interaction 
is included in QRPA changing the coupling constant / from /=0.0, 0.5, 1.0, 1.5 
and 1.7 in Eq. (2)- The solid lines show the unperturbed strength without RPA 
correlations. The QRPA strength is smoothed out by a Lorentzian function with 
a width T = 1 MeV 

The density dependent contact pairing interactions are 
adopted for both T=l and T=0 channels, 

Vr=i(ri,r 2 ) = Vo^^fl - — Wi - r 2 ), (1) 
2 V p I 

1+Po-/ P(r) \ 
Vr=o(ri,r 2 ) = fV —^ 1 - ^ Wr, - r 2 ), (2) 
2 V p I 

where r = (ri + r 2 )/2 and po is the saturation density taken 
to be po=0.16 fm -3 . The strength of T=l pairing is deter- 
mined to reproduce the average odd-even mass staggering gap 
A (3) = 12/A 1/3 MeV in N=28 region to be V =-550MeV-fm 3 . 
We solve firstly the HFB equations in coordinate space with 
a large radial mesh extending up to 20fm with a step of 0.05fm. 
In QRPA, we take into account all the states up to 180MeV ex- 
citation energy which have the occupation probabilities more 
than v 2 = 10~ 6 . The T=0 pairing is changed from null pair- 
ing to strong pairing limit taking the factor / from /=0.0, 0.5, 
1.0, 1.5 and 1.7. The f=1.7 is the maximum value to avoid the 
spin-isospin instability of the ground state of 55 Cu. 

The GT transition involves both the spin and isospin opera- 
tors as 

0{GT) =Y u o-{i)Ui) (3) 

i 

for spin-flip and isospin-flip excitations. The calculated GT 
states in 56 Ni for f_ channel are shown in Fig[T]where the inter- 
action SGII+Tel 1I24I.I25I1 is used for HF mean field and particle- 
hole interactions. The unperturbed strengths shown by solid 
straight lines have two groups at E A =9 MeV and 17 MeV. The 
QRPA results without T=0 pairing shown by the black curve 
have the main peak around 20.4 MeV together with a small 
peak at 18.0 MeV. When the T=0 pairing strength is introduced, 
the global GT strength is shifted to lower energy. At the same 
time, the GT strength of lower peak is increased and absorbs 



Table 1: Amplitudes of main (np) particle-hole and particle-particle type con- 
figurations of GT states in 56 Ni. The QRPA calculations are performed without 
and with the T=0 pairing interaction in the cases of / = and / = 1.5, respec- 
tively. The Skyrme interaction T21 is used for HF and p-h matrix calculations. 
The abbreviations B and C correspond to the GT reduced matrix element and 
the normalization factor defined in Eqs. |4j and {5), respectively. The excitation 
energy Ex is given in unit of MeV . 



'Ni / = 



E, 


B(GT) 




B C 


17.5 


1.41 


2p 3/2 (0.20) 2p 3/2 (0.21) 
lf 7/2 (0.69) lf 5/2 (0.09) 
lf 7/2 (0.69) lf 7/2 (0.67) 


0.358 0.127 
-0.229 0.006 
1.341 0.802 


18.3 


1.49 


2 Pl/2 (0.10) 2p 3/2 (0.21) 
2p 3/2 (0.20) 2p 1/2 (0.11) 


0.260 0.153 
0.846 0.740 


21.3 


7.88 


lf 7/2 (0.69) lf 5/2 (0.09) 


2.48 0.742 



S_(GT)=18.28 



'Ni /=1.5 



E, 


B(GT) 




B C 


16.6 


4.82 


2p 3/2 (0.20) 2p 1/2 (0.11) 
2p 3/2 (0.20) 2p 3/2 (0.21) 
lf 7/2 (0.69) lf 5/2 (0.09) 
lf 7/2 (0.69) lf 7/2 (0.67) 


-0.203 0.049 
-0.682 0.491 
-0.237 0.007 
-0.790 0.339 


20.5 


4.10 


lf 7/2 (0.69) lf 5/2 (0.09) 


-1.900 0.437 



S_(GT)=14.75 



the major strength when the T=0 pairing is stronger than T=l 
(/>!)■ 

The HFB+QRPA results for ?_ channel are shown in Fig. [2] 
for 48 Cr, 56 Ni and 64 Ge. The interaction T21 is used for HF 
mean field and particle-hole interactions. When the T=0 pairing 
strength is introduced, the global GT strength is shifted to lower 
energy. At the same time, the increase of lower energy peak 
is always seen in the three nuclei as the T=0 pairing becomes 
stronger. As a systematic trend, the lighter the mass is the larger 
the lower energy GT peak is. In the two nuclei 56 Ni and 54 Ge 
, the double peak structure can be seen even in the strong T=0 
pairing limit at / = 1.7. On the contrary, in 48 Cr, the higher 
energy peak is almost disappeared in the case of / = 1.7. We 
found this trend also in 44 Ti in which the higher energy peak is 
very small and is not seen in the case of / > 1.5. 

The experimental data of 56 Ni show two peaks at 18.7 and 
20.8 MeV H20TI . as is shown in the middle panel of Fig. 2. The 
calculated results for 56 Ni with T=0 pairing of / = 1.0 and / = 
1 .5 give almost the identical peak energy to the observed one 
for the higher excitation peak. On the other hand, the calculated 
lower peak is 2 MeV lower than the experimental one in energy. 

Main configurations of dominant GT states in 56 Ni are tab- 
ulated in Table Q] without and with the T=0 pairing interaction 
(/ = 0.0 and / = 1.5), respectively. In Table [1] the value B 
corresponds to the GT reduced matrix element 

B = (Xu K v v - Yu v v n ){n\\6{GT)\\v) (4) 

with the operator Eq. (0, and the value 

C = X 2 - Y 2 , (5) 
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is the normalization factor where X and Y are QRPA ampli- 
tudes. The dominant configurations are only listed in Table 
1 . Without the T=0 pairing, the main configuration of higher 
excitation is a particle-hole type (lf^ 2 ,lf^ 2 ) excitation refer- 
ring to the configuration of fully occupied If 7/2 proton and neu- 
tron orbits, while that of lower energy has particle-particle type 
(lf^ 2 ,l/^ 2 ), (2pL 2 ,2pj, 2 ) and (2p^ 2 ,2p^ 2 ) excitations. These 
particle-particle type excitations become possible only when 
T=l pairing interaction acts on the ground state. The occu- 
pation probabilities of neutron single particle states are ob- 
tained by integrating the lower components of HFB wave func- 
tions as v 2 =0.69, 0.20 , 0.10 and 0.08 for lf 7/2 , 2p 3/2 , 2p i/2 
and IT5/2 orbits, respectively. The proton states have almost 
the same occupation probabilities. With a strong T=0 pairing 
(f-1.5), the main configurations of the low energy peak be- 
come a coherent superposition of two quasi-particle (QP) exci- 
tations (2p* /2 ,2p* /2 ), (2p^ /2 ,2pf /2 ), (lf? /2 ,lfJ /2 ) and (lf^ /2 ,l# /2 ) 
as shown in Table [TJ These coherent phases are inherent to 
the T=0 pairing interaction and increases the GT strength. On 
the other hand, the higher energy peak is still dominated by 
the (1£^ 2 , lfi^) configuration, but decreases its strength signif- 
icantly. 

The GT sum rule values S_(GT) are given in Table[TJwithout 
and with the T=0 pairing (/ = 0.0 and / = 1.5). The S_(GT) 
value is decreased when the T=0 pairing is included [21]. This 
is due to the fact that the high excited GT strength is more de- 
ceased than the increase of the lower excited GT strength since 
the main configuration of GT strength (l£™lf^ 2 ) lS hindered. 
The decrease of GT strength is always the case with T=0 pair- 
ing also in the other N=Z nuclei studied in the present paper. 
The decrease of the S_(GT) value is about 20% and going up to 
30% in the cases of / = 1 .5 and / = 1 .7, respectively. We have 
checked the increase of GT strength in the lower peak in 56 Ni 
by using other Skyrme interactions SGII+Te 1 1124 12511 . T43 and 
T32 [26]. We found that these features of GT states are robust 
with some minor differences in the excitation energy of two 
main peaks. 

The enhancement of lower energy GT strength by the (T=0, 
S=l) pairing can be induced by the transition between four QP 
states such as 

<{(/)s=0(P3/2P3/2)^)^^^ 

= ((P3/2P3/2) S r =?l|0(Gr)||^^ /2 ^^2)^) (6) 

where the two quasi-particle configuration (j 2 ) T=l ' s=0 (for ex- 
amples, j - lfy 2 or J ~ 2py 2 ) acts as a spectator. 
The double bar in Eq. (0 is the symbol of doubly re- 
duced matrix element both in spin and isospin space. Other 
QP configurations (l/ 7 " /2 , lf 7 ) 2 ) T=0 - s =\ (l/ 7 v /2 , 1% /2 ) T=0 ' S=1 and 
(2py 2 Zpy 2 ) T=0 ' S=i are also possible configurations in Eq. © 
equivalent to (2p* 2 2py 2 ) T S=1 . It should be noticed that 
the T=l pairing works between the two particle configurations 
|(Z;) 2 : (L = S = 0)J = 0,T = 1) with j = 7/2, 5/2, 3/2 and 
1 /2 in pf shell model space. On the other hand, the T=0 pair- 
ing acts on the configurations \(lj,l'f) : (L = 0,5 = l)J = 
1,7 = 0) with not only (Z = j = /), but also with (/ = I', 
j = f ±1). This is the reason why the T=0 pairing makes 
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Figure 2: (Color online) Gamow-Teller strength in 48 Cr, 56 Ni and 64 Ge by 
HFB+QRPA with the Skyrme interaction T21. The quenching factor 0.74 is 
adopted for the GT transition operator O(GT). The T=0 pairing interaction is 
included in QRPA calculations by changing the coupling constant from /=0.0, 
0.5, 1.0, 1.5 and 1.7 in Eq. <2) The experimental data of 56 Ni are taken from 
Rot. i20h. 



strong couplings among (T=0,S=1) pairs with (2py 2 2p^ 2 ) and 

( l fi//2 l f5p) aS Wel1 aS ( 2 ^3/2 2 ^3/2) and ( l fll2 l fll2>> enhancin g 

the lower energy GT peak in the case of / = 1 .5 in TableQ] 

Main configurations of GT states in 48 Cr and 64 Ge are listed 
in Table [2] calculated with a strong T=0 pairing (/ = 1.5). Be- 
cause of partially occupied 2p3/ 2 , lf7/2 and ld3/ 2 orbits in 48 Cr, 
the configurations (2p^ /2 2p5 /2 ), (lfi; /2 lf* 2 ) and (ld^ld^) 
give large contributions to the lower energy GT strength at 
E t =17.4 MeV having coherent phases. Although the high en- 
ergy state at E x =19.4 MeV obtains large contributions from 
(2p^ 2 2p^ 2 ) an d (l r 7/ 2 l r 5/ 2 )' ^ i s cancelled by other config- 
urations (ld^ld^j) an d ( ir 7/ 2 ir 7/ 2 )- This cancellation de- 
creases substantially the B(GT) value compared with that of 
the state at E A = 17.4 MeV. The situation is rather different in 
64 Ge since neutron and proton lf?/ 2 orbits are almost fully oc- 
cupied. Thus the (lf^ 2 lf^ 2 ) configuration does not appear as 
a dominant configuration in the strong B(GT) states in Table 
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II. On the other hand, the 2pi/2, 2p3/2 orbits contribute signifi- 
cantly to the lower energy state at E x = 16.1 MeV in 64 Ge. The 
B(GT) value of this state in 64 Ge is only a half of the state at 
E A = 17.4 MeV in 48 Cr because of negligible contribution from 
the (lf^, 2 lf^, 2 ) configuration even though the T=0 pairing is 
strong. The high energy state at E A =21.2 MeV in 64 Ge is dom- 
inated by the (lC^li^) configuration, the same as the high 
energy state in 56 Ni. The sum rule values S_(GT) of 48 Cr and 
64 Ge are also given in Table|2]in the case of strong T=0 pairing. 
In general, the sum rule value is larger for the heavier N=Z nu- 
clei since available configuration space is larger. The increase 
of S_(GT) is 25% and 47% in 56 Ni and 64 Ge, respectively, in 
comparison with that of 48 Cr. 



Table 2: Same as Tabled but for 48 Cr and 64 Ge with the T=0 pairing interaction 
/ = 1.5. 



s Cr 



/= 1.5 





B(GT) 


v(vj) 




B 


C 


17.4 


5.68 


2p 3 / 2 (0.12) 


2p 3 /2(0.12) 


0.186 


0.062 






ld 3/2 (0.85) 


ld 3 / 2 (0.84) 


0.251 


0.204 






lf 7 / 2 (0.33) 


lf 5 / 2 (0.07) 


0.268 


0.021 






lf 7 / 2 (0.33) 


lfv/2(0.32) 


1.04 


0.558 


19.4 


1.19 


2p 3/2 (0.12) 


2pi /2 (0.07) 


0.215 


0.081 






2p 3 / 2 (0.12) 


2p 3 /2(0.12) 


0.559 


0.461 






ld 3 / 2 (0.85) 


ld 3/2 (0.84) 


-0.217 


0.142 






lf 7 / 2 (0.33) 


lf 5 / 2 (0.07) 


0.383 


0.038 






lfv/2(0.33) 


lfv/2(0.32) 


-0.384 


0.054 



S_(GT)=11.77 



64 



Ge 



/ = 1.5 



E, 


B(GT) 


v(vj) n(wl) 


B C 


16.1 


2.15 


2p 3/2 (0.48) 2p 1/2 (0.21) 


-1.316 0.889 


21.2 


5.21 


lf 5/2 (0.14) lf 7/2 (0.90) 
lf 7/2 (0.92) lf 5/2 (0.15) 


-0.187 0.353 
2.392 0.561 



S_(GT)=17.26 



The peak ratio between GT strength of the low energy region 
(below E v =20 MeV) and the high energy region (above E A -= 20 
MeV) is plotted in Fig. [3] as a function of mass number A. In 
the cases of weak T=0 pairing (/ < 1 .0), the high energy peak 
Bhigh(GT) is larger than the low energy one Bi ow (GT). With 
stronger T=0 pairing (/ > 1 .0), the Bi 0W (GT) value becomes 
larger and dominates the GT strength distribution, especially in 
nuclei A < 56. The low energy peak may appear as a single gi- 
ant resonance in the strong T=0 pairing limit. The empirical ra- 
tio of B] ow (GT) to Bhigh(GT) in 56 Ni is obtained from observed 



B(GT) strength distributions in Ref. 12011 . The data point is con- 
sistent with the calculated ratio with the T=0 pairing at / ~ 1 .5. 
The shell model calculations were performed for 56 Ni by us 



l by us- 
ing effective interactions KB3G and GXPF1A in ref. 
and a particle-vibration coupling model is applied to calculate 
the GT strength in 56 Cu in ref. 12811 . The shell model interaction 
GXPF1A gives a two-peak structure for GT strength distribu- 
tions. The competition between T=0 pairing and the spin-orbit 
interaction may be a key issue of these results 12711 . It is also 
shown in ref. [28] that the effect of particle-vibration coupling 
gives a sizable effect on GT strength distribution in 56 Cu based 
on the HF+RPA model. It might be an interesting and chal- 
lenging task to take into account the effect beyond HFB+QRPA 
such as the particle-vibration coupling effect in the microscopic 
model in future. 

The strength of the T=0 pairing can be determined by using 
the p-n scattering length H29I1 in the same way as for the T=l 
pairing ||30[ 31 ] . The extracted value is slightly larger than the 
T=l pairing strength and, makes a bound deuteron state with 
a binding energy of 1.41 MeV, which is somewhat smaller than 
the experimental value 2.2MeV. The effective density depen- 
dent pairing forces of zero range are obtained in T=0 and 1 
channels in symmetric nuclear matter from Paris forces in Ref. 
Ill4ll . The obtained T=0 interaction is much stronger than the 
T=l interaction. In Ref. 113211 . a phenomenological T=0 pair- 
ing strength was determined to be 1.36 times stronger than the 
corresponding T=l pairing by fitting the ground state spins of 
several N=Z nuclei in the p shell region. The large value of the 
T=0 pairing can be considered to mimic the effect of tensor cor- 
relations which contributes to the binding energy of deuteron 
significantly. At present, it is still an open question how much 
is the phenomenological T=0 pairing interaction, which can be 
used for HFB and QRPA calculations in medium-heavy and 
also in heavy nuclei. As we have pointed out in this letter, the 
empirical data of GT strength in N=Z nuclei will provide im- 
portant information to disentangle the T=0 pairing in medium- 
heavy and heavy mass nuclei. 

In summary, we study GT states in several N=Z nuclei 
in medium mass region by using HFB+QRPA model taking 
into account T=0 pairing interaction. We found in the con- 
text of SU(4) supermultiplets that the low excitation energy 
GT strength is largely enhanced by T=0 pairing interaction 
which is working cooperatively with the T=l pairing corre- 
lations in the ground state. Especially in the lighter nuclei 
N=Z=24~28, the Bi ow (GT) becomes comparable or even larger 
than the Bh; g h(GT) values in the strong T=0 pairing (/ > 1.0). 
This enhancement of the lower energy GT strength is attributed 
to the coherent excitation of several pf shell configurations and 
can be considered as an implementation of the super-allowed 
GT transition in the SU(4) supermultiplets in the spin-isospin 
space. It is found that the available data in 56 Ni is consistent 
with the calculated ratio Bi ow (GT)/Bhi g h(GT) with the T=0 pair- 
ing strength / ~ 1.5. Further additional information of empir- 
ical B(GT) is desperately desired to confirm the role of T=0 
pairing in the GT states in N=Z nuclei. 

We are grateful to Toshio Suzuki, M. Honma, M. Ichimura, 
K. Yako and T. Wakasa for valuable discussions. This work 
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Figure 3: (Color online) Ratio between Gamow-Teller peaks in the lower en- 
ergy and the higher energy regions as a function of the T=0 strength parameter, 
/, calculated with T21 interaction. The lower energy and higher energy re- 
gion is divided at E v =20MeV. The experimental ratio in 56 Ni is obtained from 
observed B(GT) strengths in the lower and higher energy regions in Ref. [20] . 
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